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When analyzing the particle-like excitations in quantum field theory it is natural to regard 
the field mode corresponding to the particle momentum as an open quantum system, together 
with the opposite momentum mode. Provided that the state of the field is stationary, 
homogeneous and isotropic, this scalar two-mode system can be equivalcntly represented in 
terms of a pair of quantum Brownian oscillators under a Gaussian approximation. In other 
words, the two- mode system behaves as if it were interacting linearly with some effective 
environment. In this paper we build the details of the effective linear coupling and the 
effective environment, and argue that this quantum Brownian representation provides a 
simple, universal and non-perturbative characterization of any single particle-like excitation. 
As immediate applications of the equivalence, we reanalyse the interpretation of the self- 
energy in terms of decay rates in a general background state, and present the master equation 
for the field mode corresponding to the particle momentum. 



I. INTRODUCTION AND MOTIVATION 

Quantum fields can be regarded from the viewpoint of open quantum systems [TJ [2j |3j 0]. 
The degrees of freedom of the field which are relevant for the physical problem in question 
constitute the reduced subsystem and the rest form the environment. If there are several 
fields in interaction, and the object of interest is a particular field, it is natural to trace 
over the other environment fields. For instance, in electrodynamics one can study the so- 
called Euler-Heisenberg effective action for the photons [5] , considering the electrons as the 
environment, or the complementary case, in which the electrons are taken as the system 
of interest and the photons are integrated out (6j El El [9] . Similarly, in stochastic gravity 

the system of interest is the gravitational 
field and the matter fields are integrated out. In many other circumstances it is natural 
to consider as the reduced system the modes of the quantum field which are below some 
ultraviolet cutoff, with the ultraviolet modes constituting the environment. This approach 
has been used, for instance, studying bubble nucleation [221 [23] , analyzing decoherence in 
field theory [2H 123 [2S] and in inflationary cosmology [23 EH1 1221 ED] • As these references 
illustrate, the open quantum system point of view has often provided new tools and insights 
to different field theory problems, specially when dealing with states different than the 
Minkowski vacuum. 

The propagation of a particle-like excitation over a given field background is another 
situation in which there is a natural system-environment separation: the field mode corre- 
sponding to the particle momentum is the object of interest, and the rest of modes of the 
field, along with any other field in interaction, form the environment. In the Minkowski 
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vacuum, the Kallen-Lehmann representation of the propagator l31"| 132] [33] provides a 
complete description of single particle excitations from the field theory perspective, but 
there was no equivalent analysis in thermal backgrounds [Ml [351 133 EH EE] , oi in curved 
spacetimes [391 HO) HJJ [32] . The open quantum system approach proves a useful tool when 
reanalising the propagation of particle-like excitations in non-vacuum or non- Minkowski 
backgrounds [HI HI] . 

Therefore, in this paper we elaborate on the open quantum description of the field mode 
corresponding to the particle momentum, with the goal of developing useful techniques for 
the analysis of particle and quasiparticle excitations. We concentrate on the basic methods 
and tools, and give example applications; the general analysis of the particle-like excitations 
in non- vacuum or non-flat backgrounds is left for separate publications |43[ HI] . 

We would like to consider the field mode as the system and the remaining modes as the 
environment. However, it proves difficult to implement this system-environment division 
directly, since the mode-decomposed field (ft p is not real, but is a complex quantity obeying 
the contraint (ftp = <ft*_ p and acting on both the Hilbert space sector with momentum p and 
the sector with momentum —p. Instead of focusing on a single mode, given that the field 
naturally links modes with opposite momentum, we shall choose as the system of interest 
any two modes with a given opposite momentum, and as the environment the remaining 
modes of the field, as well as the modes of any other field in interaction. Namely, given 
a single scalar field (ft, the system degrees of freedom are the two modes (ftp and <ft- p , and 
the other modes (ftq, with q ^ ±p, form the environment. A similar system-environment 
division has been used in Refs. [45 [ HH1 EZ] m an inflationary context. 

The main goal of this paper is to show that this two-mode system behaves, under certain 
assumptions, as if it was an open quantum system interacting linearly with some effective 
environment, even if the original field interaction is non-linear. The paradigm of linear open 
quantum system is the quantum Brownian motion (QBM) model, whose system of interest 
is a non-relativistic particle interacting linearly with an infinite bath of harmonic oscillators. 
This model has had many applications in different contexts, among which one may mention 
the quantum to classical transition [H] H9] , the escape from a potential well [501 ED E21 [53] , 
the Unruh effect [Ml 55J or quantum optics [3] [56] . In an influential paper Caldeira and 
Leggett [57] applied the influence functional model of Feynman and Vernon |58[ 159] to the 
QBM model. The QBM model can be generalised to encompass the general class of linearly 
interacting open quantum systems |6LH I61|. 

The results of this paper can bee seen as an explicit implementation of the statement 
by Hu and Matacz [55] that the motion of a Brownian particle can be used to depict the 
behaviour of a single quantum field mode. 

We will allow for stationary and isotropic, but otherwise arbitrary, states for the field 
(a Gaussian approximation will be also assumed later on). Field theory with arbitrary 
field states can be studied within the closed time path (CTP), or in-in, method, which was 
originally proposed by Schwinger [62] and Keldysh [63]. The most characteristic feature 
of the CTP method, in contrast to the conventional in-out method, is the doubling of the 
number of degrees of freedom. 

The paper is organised as follows. In Sect. [IT] the relevant system-environment separation 
is discussed. In Sect. Ill, the central section of this paper, we present and analyse the QBM 
representation for the field modes, given the structure of the generating functional and the 
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two-point propagators. In Sect. IV we illustrate the utility of the QBM representation 



with two example applications. First, we rederive the interpretation of the imaginary part 
of the self-energy. Second, we build and analyze the master equation for the field modes 
corresponding to the particle momentum. The main body of the paper ends with Sect. |VJ 
where we summarise and discuss the main results. In order not to break the continuity of 
the exposition, some details of the analysis are left for the appendices, which also provide 
background reference material which fixes the notation and makes the paper relatively self- 
contained. In appendix |A| we briefly present the CTP approach to field theory and apply it 
to the analysis of the structure of the two-point propagators, and in appendix [B] we introduce 
the theory of linear open quantum systems, focusing also in the structure of the two-point 
propagators. The analysis of the propagators done in the appendices is important for the 
discussion in Sect. IIIII 

Throughout the paper we work with a system of natural units with h = c = 1, denote 
quantum mechanical operators with a hat, and use a volume-dependent normalisation in 
the definition of the field modes [see Eq. ||5) below]. The same symbol will be used for a 
quantity and its Fourier transform whenever there is no danger of confusion. 



II. FIELD MODES REGARDED AS OPEN QUANTUM SYSTEMS 

Let us now present the relevant system-environment separation. For concreteness, we 
consider a self-interacting field theory model consisting of a single scalar field eft, although 
results can be straightforwardly extended to any number of fields. The field <fi can be 
decomposed in modes according to 

<Ap=-^=| d 3 xe-**0(x), (1) 

where V is the volume of the space, a formally infinite quantity. The factor V~ l l 2 in the 
definition of (f> p is chosen so that the propagators verify: 

G+(M';p) = J d 3 x e-^Ol^x^i'^lO) = <O|0 p (^_ p (t')|O). 

As stated in the introduction, given a particular momentum p / 0, the system is com- 
posed by the two modes (ft p and <^>_ p , and the environment is composed by the other modes 
of the field, c/> q , with q 7^ ±p. Should there be other fields in interaction of any arbitrary 
spin, the modes of these additional fields would also form part of the environment. 

The Hilbert space can be decomposed as TL = TC sys 'Si'H env , where in turn 7i sys = TL\ i ®'H-\ i . 
Notice that this separation does not correspond to the Fock space decomposition. The entire 
system is in a state p; the state of the reduced system is p s = Tr env p, and the state of the 
environment is p e = Tr sys p. Generally speaking, the state for the entire system is not a 
factorised product state (i.e., p ^ p s ® p e ). 

The action can be decomposed as S = S sys + S^ount + S'env + SW, where S sys is the 
renormalized system action, 

= dt (j> p c/>-p - Ep0p^_pj , (2a) 
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S^ount is the appropriate counterterm action, 



S, 



count 



J dt {(Z p - 1)4> P 4>- P - [Z p (p 2 + m 2 ) - £?]M-p} > ( 2b ) 



(with mo being the bare mass of the field), and S env and S{ n t are the environment and 
interaction actions, respectively, which depend on the particular field theory model. Notice 

1 /2 

that we have allowed for an arbitrary rescaling of the field (j) — > <ft/Z p and for an arbitrary 
frequency of the two-mode system E p [which needs not be necessarily of the form (p 2 + 
m 2 ) 1 ' 2 ]. 

Let us draw our attention on the field rescaling and the frequency renormalisation. Since 
it is always possible to freely move finite terms from the system to the counterterm action and 
vice versa, both the field rescaling and the frequency renormalisation should be taken into 
account even if no infinities appeared in the perturbative calculations. A physical criterion 
needs to be chosen in order to fix the values of these two parameters. In the vacuum, such 
criterion is provided by the on-shell renormalisation scheme; in non-vacuum situations, it is 



investigated in Ref. [13], and will be briefly discussed in Sect. IV Notice that the form of 
Z p and E p is not necessarily dictated by the Lorentz symmetry: even if the countertems 
which remove the infinities from the vacuum theory also remove the infinities in general field 
states, there can be finite Lorentz-breaking contributions. Anyway for the purposes of most 
of this paper the values of Z p and E p are not be relevant and will be left unspecified (some 



comments will be made in Sect. IV B though). 

The system variables <j) p and <j)- p are complex quantities verifying <jf_ p = P . We can 
construct real degrees of freedom by introducing the following change of variables: 

1 —i 

0E = ^Op + ^-p), = ;^=(0p-0-p)- (3) 

In terms of the real variables, the system action can be reexpressed as 

S sys = i J dt (fa - E%<\>\ + - £ 2 <M) . (4) 

We could alternatively have directly obtained these two real degrees of freedom by work- 
ing with the sine and cosine Fourier transform [55]. However we prefer to work with the 
exponential Fourier transform to make manifest the momentum conservation properties. 

Most information on the reduced quantum system can be extracted from the set of 
correlation functions, or equivalently from the CTP generating functional for the reduced 
system, which can be written as 



Z\ja,a] = exp 



1 

~ 2! 
1 



4! 



dtdt' m^G^t') 

(5) 

dt dt' dt" dt'" £(t)j^')^(i")^(t"')G2^(i, t', t", t'") + 



This expression is somewhat cumbersome and needs some clarification. Latin indices 
a,b,c. . . are CTP indices and take the values 1 and 2, indicating respectively the forward 
and backwards time branches characteristic of the CTP formalism (see appendix |A[) . Greek 
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indices a, /3, 7, . . . take the two values +p and — p, and make reference to the two field 
modes 4> p and <f>- p . An Einstein summation convention is used both for Latin and Greek 
indices. The propagator G^(t,t') is the 2-point propagator connecting CTP indices a and 
b, whose external legs correspond to particles with momenta a and (3. When the state is 
translation-invariant and isotropic, momentum conservation imposes: 

Git P)( " P) (M') = Gi/ K+P) (M') = G ab (t,t'; P ), 
Git PK+P \t,t') = G { - b p)i ^(t,t') = 0. 

In turn, G^2?^ 7(5 \t,t' ,t" ,t"') is the connected part of the four-point correlation function 
having external legs with momenta a, (3, 7 and 5. For translation-invariant states momen- 
tum conservation implies that only when momentum is balanced (i.e., two incoming and 
two outgoing external legs) the correlation function is non- vanishing. Terms with a higher 
number of external legs behave similarly. 

The open quantum system is non-linear, and a systematic treatment of the generating 
functional can be done by using the tools of non-linear open quantum systems. We follow a 
different path in the next section. 



III. QUANTUM BROWNIAN MOTION ANALOGY 



Let us reconsider the generating functional, Eq. (|5j . It depends on the n-point correlation 
functions, with n being arbitrarily large. However, in many situations one is interested in 
properties which only depend on the two-point correlation functions. In other cases one 
is doing a perturbative expansion of the generating functional, and connected higher order 
correlation functions are usually also of higher order in the expansion parameter. Finally 
there are situations in which one only has access to the two-point correlation functions, and 
expects (or simply hopes) that connected higher order correlation functions are subdominant. 
In any of this situations one can be tempted to approximate the generating functional by 
the following Gaussian expression: 



Z[ja,a] ~ exp 



1 

2! 



dtdt'f a (t)j b p(t)G a J(t,t>) 



(7) 



This Gaussian approximation can be controlled in the framework or the large- N expansion 
(where TV" is the number of scalar fields) [Ml ESI ES] • We should emphasise that it does not 
necessarily imply any perturbative expansion in the coupling parameter, nor any free field 
approximation. 

The Gaussian approximation of the generating functional can be expanded to 



Z[ja,p,ja-p] = exp 



dtdt' 



jp(t)j b - p (t')G ab (t, t') + j°(t')jt p (t)G ab (t, t') 



(8) 



or equivalently, in terms of the real variables 0s and </>a and its corresponding classical 
sources, 



Z\ja,T.,ja,A\ = ex P 



dtdt' 



m)iUt')G ab (t, t') + j a A (t)j b A (t')G ab (t, t') 



O) 
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Original system 



Equivalent linear QBM 



System 
System d.o.f. 
Environment 
Env. d.o.f. 
Frequency 
Coupling 
2-point function 



two field modes 

0p, </>- P 
other modes 3-d field 
</> q , q^±p 



A y- 

3! ^qq' 



J P 

*>±p0q0 q " 

G a6 (i,t';p) 



two identical oscillators 
2 copies of q 
2 1-d field 

n(=) 



G ab (t,t') (=) 



TABLE I: Detail of the equivalent linear QBM system for a A(/> 4 quantum field theory. The A0 4 
model has been chosen for concreteness, but the correspondence would be analogous for any other 
field theory model. The symbol (=) indicates that the original and equivalent quantities are indeed 
identical despite the name change. 



or 



Z\ja,vJa,A\ = Z[j a ^]Z[j aA ], Z[j a ] = exp 



dtdt'j a (t)j\t')G ab (t,t') 



(10) 



Eq. (10) shows that for translation-invariant states within the Gaussian approximation the 



reduced two-mode state effectively behaves as two decoupled quantum mechanical degrees 
of freedom. 

The Gaussian approximation implies that the expression of the generating functional of 
a reduced two-mode in terms of the two-point propagators coincides with that of a QBM 



model, provided the system is isotropic and translation-invariant: compare Eq. (B24) with 
Eq. (10). Moreover, by comparing Eqs. (B23) with Eqs. (A19), we also realise that the 



structure of the two-point propagators is identical in both cases. Notice that this latter fact 
is independent of the Gaussian approximation. 

Therefore, we conclude that, assuming homogeneity, isotropy and stationarity, there is 
an equivalent QBM system for every scalar two-mode pair treated under the Gaussian ap- 
proximation. In other words, within the Gaussian approximation, every two-mode of a given 
quantum field theory can be described in terms of a pair of quantum Brownian particles 
interacting linearly with some effective environment. We must stress that, similarly as the 
linear interaction does not coincide with the real coupling, the effective environment does not 
coincide with the real environment. The precise details of the equivalence are summarised 
in table [i] for the particular case of the A<^> 4 theory. We will see below that both the strength 
of the linear coupling and the state of the equivalent environment depend on the details of 
the original environment. 

Let us investigate on this correspondence. On the one hand, as is is shown in appendix 
|B"| the effect of the environment in the QBM system is fully encoded in two kernels, the 
dissipation kernel D(t,t') — or its closely related counterpart H(t,t'), see Eq. (Bllb) — and 
the noise kernel N(t, t'). The dynamics of the quantum Brownian particle can be determined 
once the frequency and the noise and dissipation kernels are known. On the other hand, as 
shown in appendix [Aj in a quantum field theory the two-point correlation functions are fully 
characterised by the frequency E p and the retarded and Hadamard self-energies, t') 
and EW(i, t') respectively. By comparing again Eqs. (B22) and (B23) with Eq. ( A19| ) we 
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realise that the precise analogy goes as follows: 



O — E p , 

H(t, t') = p 2 + m 2 -E p + E R (t, t'; p), 
D(t,0 = -iImS R (t,t';p), 

iV(t,t') = ^( 1 )(t,t / ;p), 



(11a) 
(lib) 
(11c) 

(lid) 



where the quantities on the right hand side correspond to the original field theory system 
and the quantities on the left hand side correspond to the equivalent QBM model. The 
correspondence is valid no matter the renormalisation scheme chosen to fix Z p and R p . 

This representation provides thus a first rough interpretation for the retarded and 
Hadamard self-energies, £ R (t, t';p) and SW(t, t';p) respectively. The retarded self-energy 
corresponds to the dissipation kernel, so that it determines the dissipative properties of the 
system, and it is independent of the state of the equivalent environment (though not inde- 
pendent of the state of the original environment). The Hadamard self-energy corresponds 
to the noise kernel, and thus it is basically related to fluctuations. 

Although the description in terms of the noise and dissipation kernels is often sufficient, 
the equivalent QBM system can be alternatively described in terms of the linear coupling 
constant to the effective environment, g, the distribution of frequencies of the environment 
1{uj) (see appendix IB]), and the occupation number of the modes of the effective environment 
n(p) = Tr[p^ S ^ OpCLp]. The product g 2 I(u>) : which determines the coupling strength to the 
the o;-mode of the effective environment, can be obtained from the imaginary part of the 
self-energy: 



ImS R (o;;p) = iD(u) 



(12a) 



This last equation implies that the equivalent coupling depends on the state of the real 
environment, since the retarded self-energy is state-dependent in general. The occupation 
numbers of the effective environment n{p) can be reproduced from the Hadamard self-energy: 



S (1) (w;p) = -2iN(u) = -2ig 2 \uj\I(Lo) 



1 



+ n(M) 



(12b) 



The knowledge of the occupation numbers fully determine a Gaussian stationary state for 



the equivalent environment. These results follow from Eqs. (B12) and (B14). 



Note that the analogy depends only on the Gaussian approximation, so that it can be 
extended to all orders in perturbation theory. Notice also that the results derived from the 
QBM interpretation are exact for all those properties which depend only on the two-point 
correlation functions; for the properties which depend on higher order correlation functions, 
it is a correct approximation depending on the validity of the Gaussian approximation, i.e., 
depending on the relative importance of the connected parts of the correlation functions 
with respect to the disconnected parts. 
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IV. EXAMPLE APPLICATIONS 
A. Interpretation of the self-energy 

As an first example application, let us make use of the QBM correspondence to analyse 
the physical significance of the self-energy in general backgrounds. Our findings will coincide 
with the result by Weldon [67] (which nowadays is a textbook result [37]) as far as the 
imaginary part of the self-energy is concerned. However, while Weldon's original analysis 
was only valid to first order in perturbation theory, our technique is be valid to all orders. 
Moreover, our analysis will not be tied to any field theory model. Additionally we will also 
obtain an interpretation for the other components of the self-energy. In this case the QBM 
analogy is exact since no four-point correlation functions are involved. 

To start we consider the probability that an excitation of energy uj decays into the 
one-dimensional environment, in the equivalent QBM system. The probability T_ that an 
excitation of the Browian particle with positive energy uj decays into the environment is 
given by (see e.g. Ref. [57] ) 

1 r H k 

T-H = - J ^^8{u>-\k\)\M\ 2 [l + n(\k\)l (13) 

where M is the amplitude of the transition and n(| A;|) is the occupation number of the 
environment states with energy \k\. The factor 1 + n(|fc|), which is due to the Bose-Einstein 
statistics, enhances the decay probability to those states which are already occupied. Since 
the equivalent QBM system is linear, the squared decay amplitude is simply given to first 
order in the linear coupling constant g by 

\M\ 2 = 9 2 1(uj)uj 2 , (14) 

where I(uj) is the distribution of frequencies of the effective environment. The factor uj 2 
is a consequence of the derivative coupling in the QBM model. The decay probability is 
therefore: 

T_(uj) = ±g 2 l(uj)[l + n(uj)}. (15) 

Likewise, the probability that an excitation of positive energy uj is created spontaneously 
from the environment is given by 

1 f dk 1 

T+H = ^ J - \k\)\M\ 2 n(\k\) = -g 2 l(uj)n(uj). (16) 

In the original system, T_ can be interpreted as the probability that a (possibly off-shell) 
excitation with energy ui decays into the environment, and r + can be interpreted as the 
probability that an environment spontaneously creates an excitation with energy uj. 

Notice that the notion of decay rate in quantum mechanics is meaningful only when the 
excitations are long-lived, or, to put it differently, when the product g 2 Z{uj) is very small. 
There is an inherent uncertainty in the concept of decay rate, which can be traced to the 
time-energy uncertainty principle. Therefore it is sufficient to present results to first order 
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in g: the inherent uncertainty to the notion of decay rate in quantum mechanics is of the 
same order as the error done by neglecting higher powers of g. In any case, this does not 
mean that we are doing any perturbative expansion in the original system: the decay rate 
can be computed to any desired order in the original perturbative coupling constant {e.g., 
A in the case of the Xcp 4 theory) . 

We next analyse the self-energy components in the equivalent QBM mode. Given that 
ImSn(w) = (i/2)[E 21 (u>) — T} 2 {uj)\ (see appendix|A|), we start by analyzing E 21 (o;). Applying 
CTP Feynman rules we get 

- i£ 21 (t, if) = -{igfd t d t , J ^X(p) Tr cnv [pe6p(*)&(-ri(0] > (17) 

where (p\ v is p-mode of the environment field in the interaction picture. We have exploited 
that first order perturbation theory yields exact results for the self-energy in linear systems. 
Introducing two resolutions of the identity in the basis of eigenstates of the environment 
Hamiltonian it is a simple exercise [68] to show that the above equation can be developed 
to 



j P (t-t') 



Z 2 \t,t') = ig 2 d t d t , f°dp?MyV n \(n+l)e-Wt-«) +ne 
Jo P V L 

where p PjYl is the diagonal value of the reduced density matrix of the p-mode of the equivalent 
environment: r £i q ^p{n p \p e \m p ) = pp n S nm . Introducing the Fourier transform we get 

£ 2 V) = ig 2 u? I" dp ?M V p P ,n [(n + l)6(u -p) + n5(u; + p)} , 

and restricting to positive energies, 

S 21 (w) = -ig 2 Lul(uj)J2P"A n + 1 ) =ig 2 uZ(uj)[l + n{uj)}, uj > 0, (18) 

n 

where we recall that n{uS) = Tr [pi eS ^ at a^] = J2 n P^,n n 1S the occupation number the u- 
mode of the environment. By comparing with Eq. (15), we thus see that i proportional to 
the decay rate: 

S 21 (w) = 2iwr_(u;), uj>0. (19) 
Repeating the calculation for E 21 (w) we similarly find 

T, 12 (uj) = ig 2 ujl(u) ^2pui, n n = ig 2 u)Z(ijj)n(uj), uj > 0, (20) 

n 

having the corresponding interpretation in terms of the creation rate, 

S 12 (w) = 2iuT + (uj), uj>0. (21) 
When the energies are negative E 21 (o;) and £ 12 (u;) exchange roles. 
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The imaginary part of the retarded self-energy is therefore given by 

ImS R (u;) = ^[£ 21 H - £ 12 H] = -^ 2 u;X(u;), (22) 

and can be interpreted as the net decay rate for an excitation of energy oj — i.e., decay rate 
minus creation rate: 

ImS R ( W ) = - W [r„( W )-r + H]. (23) 



We therefore recover Weldon's result [67 . 

We can additionally get an interpretation for the Hadamard self-energy. It is given by 

EW(w) = -S 21 (w) - S 12 (w) = -ig 2 \u\l{u)[l + 2n(M)], (24) 

and is proportional to the probability of decay plus the probability of creation, 

EW(w) = -2i\u\[T-(uj) + (25) 

B. Master equation for relativistic quasiparticles 

As another example application of the correspondence, let us study the master equation 
for the second-quantized relativistic quasiparticles. A quasiparticle an elementary excitation 
propagating in some backgroung, characterized by a momentum p, an energy E p , and a 
decay rate 7 P verifying 7 p <C E p (so that the excitation is long-lived). The basic properties 
characterizing the long-time evolution of the relativistic quasiparticles can be extracted 
from the analysis of the retarded propagator corresponding to the field mode in a model- 
independent way [43j (differently to the short-time behaviour, which is model-dependent). 
In particular, the real and imaginary parts of the self-energy, when evaluated on shell, 
correspond to the physical energy and the decay rate of the quasiparticle, respectively: 

E 2 p = m 2 + p 2 + Re£ R (£ p ,p), 7p = -- i- Im£ R (£ p .p), (26) 

E p 

Using the QBM correspondence, these two equations can be equivalently written as: 

ReH(n) = 0, 7p = -^D(fi), (27) 

The value of E p (or equivalently O) is fixed by requiring that E p represents the physical 
energy of the quantum mode as would be measured by a particle detector. We refer to 
Ref. [43 (which makes use of the QBM correspondence) for further details on the field 
theory analysis of the quasiparticle excitations. 



Thus, let us study the master equation (B27) for q = {4>j:,4>a}- To this end it will prove 



useful to express the on-shell values of the noise and dissipation kernels in terms of the 
energy £1 = E p and the decay rate 7 = 7 P : 

D(Q) = -if(O) = ij, JV(fi) = ( - +nj , (28) 
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where n = n(\£l\) is the equilibrium occupation number of the mode in question. Notice that 
since the decay rate is necessarily small, the coupling constant to the effective environment 
is also small. Therefore the mode can be considered to be weakly coupled to the effective 
environment. 

In the weak coupling regime, the master equation coefficients are given by Eqs. (B28), 
and are usually divergent for short times [60]. These divergences are associated to the 
fact that the assumption of factorized initial conditions for the inital state is unphysical, 
as commented in appendix [Bj Since we are interested in the long-time behaviour of the 
quasiparticles, we study the asymptotic value of these coefficients, which is free from the 
short-time divergences. 

It is easy to show that the frequency shift vanishes: 



5tf 



f'OO I f'OO 

-2 / D(s, ti) cos tt(s - ti) ds = - Re / H(s, ti) e*^*) ds = Re H(Q) = 0. 
Jti ^ J_ 00 

(29a) 



Taking into account that the dissipation kernel is antisymmetric, the evaluation of the 
dissipative factor is also straightforward: 

-i roo -i roo 

r = - / D{s, t) sin Q(s - ti) ds = — / D(s, t) e tn( - s ~^ ds 

U h m% J -°° (29b) 

" ^ - T- 

The first diffusion factor is also easily computed, this time recalling that the noise kernel is 
symmetric: 



roo i roo 

Fh= I N(s,ti)costt(s-ti)ds = - / N(s, U) e^ 8 '^ ds 



(29c) 



2 ' V2 

The second diffusion factor requires some extra work: 



i roo i roo 

Ff = - N(s, U) sin U(s - U) ds = - Im / N(s, ti)6(s - U) e^ 8 '^ ds 

J ti J — OO 

' ^ / '^J^k = 4P I ^ iP( " )Sig ° H [l + 2n( k |)]. 



17 J 2tt Q — oj + ie £1 J 2ir — u> 

The integrand in the last equality is only signifhcantly different from zero when uj ~ 17. 
Therefore as a first approximation we may write: 

r/w_ ^ p / ^^r^ [1+Mm = ^ ReH{m + 2n) = - (29d) 

Thus the master equation can be written as: 



i-Q^p B {q,q,t) 



1 / d 2 d 2 \ 1 2 . 2 , 2 
+ -n 2 (q 2 - a 2 



i , ,, I d d \ .„ ( 1 \ , 2 



Ps{q,q',t). 
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or, equivalently, in terms of the Wigner function, as 1 

ow s = p dw s | „ 2 dw s | d P w s | o 7 (1 , 2n) Q 2 w s 

<9i ^ <9g dp ^ dp 4 <9p 2 1 



This is identical to the master equation found by Caldeira et al. [69] , corresponding to weak 
couplings. 2 

The reduced density matrix and the Wigner function are to be interpreted in terms of the 
field modes, and not directly in terms of particles or quasiparticles themselves. Qualitatively, 
for Gaussian initial conditions the solution of the master equation can be described in terms 
of the cumulants of a Gaussian distribution in the following way |48| I7L)| . The expectation 
value of the field (q) follows the trajectory of a classical underdamped harmonic oscillator 
(q + 1Q + ^ 2 Q = 0)) namely, it oscillates and slowly decays towards the origin at a rate 
7/2. However, for single quasiparticle excitations the expectation value of the field is always 
vanishing. Therefore, it is more appealing to consider the dynamics of the second order 
cumulants, and in particular the dynamics of the energy of the mode, (l/2)(g + fi 2 g). When 
perturbed by the introduction of a quasiparticle, the energy of the mode slowly decays at a 
rate 7 towards its equilibrium value 0(l/2 + n), with n being the original occupation number 
of the mode 



V. SUMMARY AND DISCUSSION 

In this paper we have explored the open quantum system viewpoint for a pair of field 
modes of opposite momentum, which are the relevant degrees of freedom for the analysis of 
particle- like excitations in field theory. The main results have been, on the one hand, showing 
that in any interacting field theory, assuming homogeneity, stationarity and Gaussianity, this 
open quantum system can be equivalently represented by two identical quantum Brownian 
particles interacting linearly with an effective environment, and, on the other hand, exploring 
the details of the equivalence, which are expressed in Eqs. (11 ) and (12) and in table [l| 



The Brownian motion equivalence is based on three simple well-known observations. 
First, the fact that the structure of the two-point correlation functions, which we reexamined 
in the appendices, is identical regardless of the nature of the interactions. A byproduct of this 
fact has been establishing a link between the self-energy and the noise and dissipation kernels, 
thus connecting the open quantum system and field theory notation and languages. Second, 
the recognition that a Gaussian truncation leads to a generating functional which coincides 
with that of a linear theory. Finally, the observation that when assuming homogeneity, 
stationarity and isotropy the two-mode pair behave as a two copies of a single degree of 
freedom. 



1 Notice that in this equation p is the canonical momentum associated to the variable q, and has nothing to 
do with the physical quasiparticle momentum p. 

2 Fleming et al. [70] found logarithmically divergent results for the master equation coefficients in the ohmic 
dissipation case, even in the asymptotic regime and in the weak coupling limit. This divergences are 
associated to the unphysical ultraviolet behaviour of the ohmic dissipation model. Our approximate 
treatment overlooks the possible divergences associated to the high energy limit of the noise and dissipation 
kernel. 
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If the system is Gaussian, stationary and homogeneous, the equivalence is exact. If the 
system is non-Gaussian, the analogy is also exact as far as two-point correlation functions 
are concerned, but only approximate for higher order correlation functions. If the system is 
non-stationary or non-homogeneous, there will be corrections to the results of the order of 
Lp, where L is the characteristic inhomogeneity time or length scale, and p is the relevant 
energy or momentum scale. Therefore the analogy is perfectly valid in non-homogeneous 
backgrounds as long as we consider modes whose characteristic wavelengths are much smaller 
than the inhomogeneity scale. 

It must be noted that for the equivalence to be useful one still needs the field theoretic 
computation of the two-point correlation functions (or, equivalently, the self-energy). The 
correspondence does not help in this calculation, but must be regarded as a tool useful for 
interpreting and analyzing the dynamics of the two-mode system. In this sense, there are 
three basic characteristics which make the QBM equivalence appealing. 

First, the equivalence is universal, in the sense that it provides the most general descrip- 
tion of the dynamics of the scalar two-mode pair within the Gaussian approximation. The 
description in terms of a linear open quantum system allows a unified description of many 
different quantum field theory systems: the details of the quantum field theory are unim- 
portant once the equivalent noise and dissipation kernels are known. Different field theory 
models can be thus classified in the same equivalence class if they lead to the same QBM 
equivalent model. 

Second, the equivalence provides a simple characterization of the interaction, given that 
an arbitrarily complicated coupling with any number of fields is reduced to a linear inter- 
action with a one- dimensional field. Linearly interacting systems have been thoroughly 
studied in the literature (see e.g. references given in the introduction and appendix [B"]) , and 
exhibit many intersting properties, among which one can cite the fact that they are exactly 
solvable. Therefore, within the regime of validity of the Gaussian approximation, all the 
methods developed in the literature for linear quantum systems can be applied to analyse 
the dynamics of a two-mode pair. 

Finally, the equivalence is non-perturbative, since it does not rely on any perturbative 
development in the coupling constant of the original field theory. It can be therefore extended 
to all orders in perturbation theory and applied to strongly interacting systems. 

The QBM correspondence, which has also been implicitly used in Refs. |71[ 172] . proves a 
useful tool when analyzing the properties of particle-like excitations in general backgrounds 
from a field theory perspective }43| . The utility of the QBM analogy is already highlighted 
by the first two immediate applications presented in this paper. 

As a first application, we reexplored the well-known result that the imaginary part of the 
retarded self-energy corresponds to the net decay rate of the particle excitations. This was 
a rather simple application which however improved the conventional textbook derivation 
by taking profit of the three properties we have remarked above: in the first place, the 
universality of the QBM analogy, given that the derivation presented in this paper was not 
tied to any specific field theory model; in the second place, its simplicity, since the calculation 
essentially reduced to doing trivial perturbative expansions in the equivalent linear system, 
and, finally, its non-perturbative character, since the derivation avoided any perturbative 
expansion in the original system. 

As a second application, we presented the relevant master equation for the dynamics 
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of the modes corresponding to the quasiparticle momentum. This was a slighly more in- 
volved application, relying on the results of Ref. [33] (and also complementing them), which 
highligted the fact that the linear open quantum system machinery can readly be exported 
to quantum field theory. The presentation and analysis of the master equation, as done 
in the paper must be understood as a first approximation to the problem, which might be 
sufficiently interesting by itself to deserve further work. 

Acknowledgments 

I am very grateful to Albert Roura for several helpful discussions and suggestions, and 
to Renaud Parentani and Enric Verdaguer for their comments on the original manuscript. 
This work is partially supported by the Research Projects MEC FPA2007-66665C02-02 and 
DURSI 2005SGR-00082. 

APPENDIX A: THE CLOSED TIME PATH METHOD AND THE TWO-POINT 

PROPAGATORS IN FIELD THEORY 

In this appendix we give a brief introduction to the closed time path (CTP) method 
(also called in-in method, in contrast to the conventional in- out method), stressing those 
aspects relevant for this paper, and apply it to the analysis of the two-point propagators. 
We address the reader to Refs. [731 Ei ES ES [77] for further details on the CTP method, 
and to Refs. [6HI Ell ESI EQl E] for further details on the structure of the two point functions. 
For the purposes of this appendix we shall consider a free or an interacting scalar field (ft, 
although results also apply for a single quantum mechanical degree of freedom. 

The path-ordered generating functional Zc\j] is defined as 

Z c [j] = Tr (pTc e <; c d*/d s xto"(*)j , (Al) 

where 4>{x) is the field operator in the Heisenberg picture, C is a certain path in the complex 
t plane, Tq means time ordering along this path and j{x) is a classical external source. 
By functional differentiation of the generating functional with respect to (ft, path-ordered 
correlation functions can be obtained: 

(A2) 

3=0 

Introducing a complete basis of eigenstates of the field operator in the Heisenberg picture, 
4>(t,x)\(ft, t) = cft(t,x)\(ft,t), as a representation of the identity, the generating functional can 
be expressed as: 

Z c \j) = J d^dV^^,^)^ (A3) 

The functional measures deft and deft' go over all field configurations of the fields at fixed 
initial time t. If the path C begins and ends at the same point ti, then the transition element 



15 



Imt 



*f Re t 



FIG. 1: Integration path in the complex-time plane used in the CTP method. The forward and 
backward lines are infinitesimally close to the real axis. 



of the evolution operator can be computed via a path integral: 

/_ _ /' l / 3 (*ii x ) =< / , '( x ) 
dWOMiW.ti) / p ¥ , e */c d */d 3 x{LM+K*)i(*)} j (A4) 

where L[(f>] is the Lagrangian density of the scalar field. 

Let us consider the time path shown in Fig. [TJ If we define yi,2(£> x ) = ¥?(£j x ) an d 
jl,z(t, x) = j(t,x) for i E Ci^, then the generating functional can be reexpressed as: 



ZbWa] = / d0dW(^ti|p|0',ti 

r^i(*f,x)=0"(x) 



x / e */d 4 a;{L[ V i]-|-ipi(a:)ji(2:)} 
Vi(*i,x)=<£(x) 
/•^ 2 (tf,x)=0"(x) 

X / 

y ¥ j 2 (t i ,x)=^'(x) 



(A5) 



In the following it will prove useful to use a condensed notation where neither the boundary 
conditions of the path integral nor the integrals over the initial and final states are explicit. 
With this simplified notation the above equation becomes 



Z[ji,j 2 ] = / VfcVfa^pW ,t) e 



! A p i I <^ x { L [fl]-L[ i P2]+V>l{x)jl(^)-'P2(x)j2(x)} 



An operator representation is also possible: 



Z\j 1 ,j 2 )=Tr[pTe 



-i J t * f dt J d 3 x fate (x) Te i /* f dt J d 3 x 4>{x)n (x) 



(A6) 



(A7) 



By functionally differentiating the generating functional, the Feynman and Dyson prop- 
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agators and the Whightman functions can be obtained: 



Gu(x,x) = G F {x,x') = TV [pT<i>(x)<j>(x')] 
G 2 i(x,x') = G + (x,x) = Tr [p^(x)^(x')] = 
Gi2(x,x') = G-(x,x') = TV [pcf)(x')(f)(x)~\ = 

G22(x,x') = Gb{x,x') = TV [f>T(f)(Q 



5 2 Z 



Sji(x)8jx(x') 
5 2 Z 



31=32=0 



Sj2(x)Sji(x') 
6 2 Z 



Sji(x)dj 2 (x') 
5 2 Z 



31=32=0 



31=32=0 



[X)0[X 



8j2(x)dj 2 (x') 



(A8a) 
(A8b) 
(A8c) 
(A8d) 



31=32=0 



These four propagators can be conveniently organised in a 2 x 2 matrix, the so-called direct 
matrix. 



, ,,_(G F (x,x') G_(x,x')\ 
Gab[X ' X) -{G + (x,x>) G»(x,x')) 



(A9) 



Lowercase roman indices may acquire the values 1 and 2 are raised and lowered with the 
"CTP metric" c a b = diag(l,— 1). Higher order correlation functions can be obtained in a 
similar way. 

We may also consider the Pauli-Jordan or commutator propagator, 



G(x,x')=Tr(p[0(x)J(x')]), 
and the Hadamard or anticonmutator function 

G (1) {x,x') = Ti(p{0{x),j ) {x')}). 



(AlOa) 



(AlOb) 



For linear systems 3 the Pauli-Jordan propagator is independent of the state and carries 
information about the system dynamics. Finally, one can also consider the retarded and 
advanced propagators 

G R {x,x') =6{x° -x'°)G{x,x') =^(x o -x /O )TV(^[0(x),0(x , )]), (Alia) 
G A {x,x') =6{x'° -x°)G{x,x') =6(x'° -x )Ty(p[4>{x),4>{x')}), (Allb) 

which also do not depend on the the state for linear systems. The retarded and advanced 
propagators and the Hadamard function can be used as an alternative basis to (A9) in the 
so-called physical or Keldysh basis. 

It is also useful to introduce the correlation functions in momentum space, which are 
defined as the Fourier transform of the spacetime correlators with respect to the difference 
variable A = x — x' keeping constant the central point X = {x + x')/2: 



G ab (u;, p;X) = j d 4 A e ^A°-i P .A G ^ x + A/2> x _ A/2) _ 



(A12) 



3 By linear systems we mean systems whose Heisemberg equations of motion are linear. These correspond 
either to non-interacting systems or to linearly coupled systems. 
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Mixed time-momentum representations of the propagator, G a ^{t, t'; p; X), can be similarly 
introduced. For homogeneous and static backgrounds the Fourier transformed propaga- 
tor does not depend on the mid point X. The canonical example of static and ho- 
mogeneous background is the thermal background, in which the state of the field is 
p = e~@ H /Tr ( e~P H ^J. Thermal field theory can be thus treated as a particular example of 
field theory over an arbitrary background. This viewpoint corresponds to the so-called real 
time approach to field theory [371 ES E2] • 

For interacting theories the self-energy can be introduced similarly to the vacuum case. 
Interaction theory mixes the two CTP branches and therefore the self-energy has matrix 
structure and is implicitly defined through the equation 

G ab (x,x') = G^(x,x>) + J d 4 yd 4 y'G^(x,y)[-i^ d (y,y')}G db (y',z), (A13) 
where G a °?(x,x') are the propagators of the free theory and G a i,(x,x') are the propagators 



of the full interacting theory. Notice that Eq. (A13) is matrix equation relating the four 
components of the self-energy with the four components of the propagator. Therefore there 
is no diagonal relation between Gn(x,x') and £ n (y, y') as in the vacuum case. The ab 
component of the self-energy can be computed, similarly to the vacuum case, as the sum 
of all one-particle irreducible diagrams with amputated external legs that begin and end 
with type a and type b vertices, respectively. CTP Feynman rules are completed with the 
prescription of adding one minus sign for every type 2 vertex. 

A particularly useful combination is the retarded self-energy, defined as Er(s, x') = 
T l n (x,x') + S 12 (x,x'). It is related to the retarded propagator through 

G R (x,x') = Gf\x,x') + f d 4 yd 4 y'G^\x, y )[-i^ R (y,y')]G R (y',z), (A14) 

This equation can be regarded as a consequence of the causality properties of the retarded 
propagator. A similar relation holds between the advanced propagator Ga(x,x') and the 
advanced self-energy Sa(x,i') = + S 21 (x,x'). Another useful combination is the 

Hadamard self-energy, which is defined as EW(x, x') = T, n (x, x')+T, 22 (x, x') [or equivalently 
as S^^Xjx') = — T, 12 (x,x') — T, 21 (x,x')] and which is related to the Hadamard propagator 
through [68] 

G W(x,x') = -i I dVV Gr(^)£ (1 W)GaO/,z')- (A15) 



All self-energy combinations can be determined from the knowledge of the Hadamard self- 
energy and the imaginary part of the retarded self-energy. This latter quantity can be 
obtained from the following cutting rule: 

ImS R (x,y) = % - [S 21 (x,y) - S 12 (x,y)] . (A16) 
So far, all expressions in this appendix refer to arbitrary background states p. For static 



and homogeneous backgrounds, Eq. (A14) can be solved for the retarded propagator by 
going to the momentum representation: 

— % 

— or + m z + p z + ^r{lo, p) — p u ie 
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We have considered that the free propagators of the mode p are those corresponding to the 
action ( 2a ) . Notice that in general the self-energy is a separate function of the energy u> and 
the 3-momentum p, and not only a function of the scalar p 2 , as in the vacuum case. The 
Hadamard function admits the following expression [which can be derived from Eq. (A15)]: 

G«(u,,p) =i|G R ( W ,p)| 2 s( 1 )(u;,p) = 2^p S v ( r P \l2^r T v < ^2' 

[-uo z + m z + p 2 + Re T, k (uj, p)\ z + [Im 2J R (u;, p)\ 2 

(A18) 

From the retarded propagator and the Hadamard function we can derive: 



, „[-o; 2 + 4 + ReS R (a;,p)]+iS( 1 )(a; ! p)/2 
Gf(w,p) = - o ft, (A19a) 

[_ W 2 + m 2 + p 2 + Re Sr(Wj p)] 2 + [Im Er(Wj p)] 2 

i [-cv 2 + El + Re E R (w, p)] + iSWfu, p)/2 

[_ w 2 + m 2 + p 2 + Re Er ( Wj p) ]2 + [Im Sr(W) p)] i 



g d(^, p) = - L , - , p , - ^\ 'L rT ' , l2 , (A19b) 



^( 1 )(cJ,p)/2+ImS R (o;,p) 
[_ w 2 + m 2 + p 2 + Re Er(Wj p)] 2 + [Im Er(Wj p) p 



G -( aJ 'P) = 7 ~ — m2 , r T ^ } v2> ^ A19C ) 



G + („, P) = iE"V.P)/2-ImE„(„,p) 

[-w 2 + m' 2 + p 2 + Re£ R (^, p)] 2 + [Im E R (ai, p)] 2 

When the field state is not exactly homogeneous, the above expressions are still correct up 
to order Lp, where L is the relevant inhomogeneity time or length scale. 



APPENDIX B: LINEAR OPEN QUANTUM SYSTEMS 



We present those aspects of the theory of linear open quantum systems relevant in this 
paper, focusing on the propagators of the Brownian particle. For a more complete presen- 
tation check Refs. [D H EH E2| • 

We shall consider a quantum Brownian motion (QBM) model: an open quantum system 
composed of a harmonic oscillator q(t), which will be the subsystem under study, linearly 
coupled to a free massless field ip(t,x), which will act as environment or reservoir. The 
action for the full system can be decomposed in the action of the harmonic oscillator, the 
action of the scalar field and the interaction term as 



Ssy a [q\ 



(it 



1 



2 

dt dx 



2 1 o2 2 

< \l q 

2 H 



Sintfe, ¥>] =9 dtdx5{x)qtf, 



(Bla) 
(Bib) 
(Blc) 



with £1 being the frequency of the harmonic oscillator and g being the coupling constant. 
The oscillator is taken to have unit mass. A counterterm action including a frequency shift 
could also be considered. 

We use a one-dimensional free field as the environment, following the treatment of 
Ref. [48 . This is equivalent to the alternative representation in which the environment 
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is modelled by a large ensemble of harmonic oscillators [57]. This equivalence can be seen 



performing a mode decomposition in the interaction term (Blc 



/ — f dp 

Smt[q,f] =VL dt—gqtp p , (B2) 



where ip p (t) is proportional to the spatial Fourier transform of the scalar field, 

<p p (t) = —j= [ dx e~ vpx cp(t,x), 



where L is the length of the real axis (formally infinite). 

The standard model, also called ohmic model, can be generalised by replacing the delta 



interaction of equation (Blc I by a function f(x). In this case the interaction term is, 

Sint[q,<p} = J dtdx f(x)gq(t)<p(t,x). (B3a) 
or equivalently in the Fourier space, 

S int [q, ( p]=VlJdt ^-f(-p) gq (t)v P (t). (B3b) 

The real even function I(p) = f(p)f(—p) is called the distribution of frequencies. The 
product g 2 Z(u>) characterises the properties of the coupling with the environment at a given 
energy u. The QBM model this way generalised encompasses the entire class of linearly 
coupled environments. 

When the system and the environment are initially uncorrelated, i.e., when the initial 
density matrix factorises — p(t\) = p s {t\) ® p e (ti), where p s (ti) and p e {ti) mean, respectively, 
the density matrix operators of the system and the environment at the initial time — , the 
evolution for the reduced density matrix can be written as 

p s (qi,q'i,t{) = J dq i dq[J{q u q^t i ]q h q[,t- l )p s {q h q[ 1 t l ), (B4) 

where the propagator J is defined in a path integral representation by 

q(t t )=q f q'(t[)=q' { 

J(q f ,q' { ,t i ;q i ,q[,t i )= J Vq J Vq' e^M-SM+SWw']) , (B5) 
5(*0=?i l'(ti)=l'i 

with Sip[q,q'] being the influence action, which is related to the the influence functional 
introduced by Feynman and Vernon |58| [59] through F[q,q'] = e lSlF ^- q ' q '\ In turn, the 
influence functional can be expressed in the following way: 

F[q,q') = JjT>( P V^p e ([<p i ]M,t l )exp[i(S[(p]-S[ V /]+ (B6) 



4 In the literature the distribution of frequencies is frequently defined as coX(u>). 
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The path integral has the boundary conditions (p(x,t{) = <fi(x), ip'(x,ti) = tp[(x), (p(x,tf) = 
ip'(x,t{) = (ff(x); there is also an implicit sum over initial and final states, <fi(x), ^[{x) and 

Considering a factorized initial state is a rather unphysical hypothesis that leads to 
surprising results in many circumstances (see for instance Ref. [60 J. The methods presented 
in this appendix can be generalized to more natural initial density matrices by the use of 
the so-called preparation functions [U [M] . However the preparation function method does 
not completely solve all the problems because it is based in a sudden change of the density 
matrix. A more physical approach involves a continuous preparation of the system [85] . In 
any case, these techniques are increasingly more involved, and we shall be only interested 
in studying the dynamics much after the typical decoherence time. In this case the system 
and environment have had enough time to interact and become entangled, and the precise 
form of the initial state becomes unimportant. 

When the initial density matrix of the environment /O e ([v?i] , [<^],ii) is Gaussian, the path 
integrals can be exactly performed and one obtains |50[ [58) 186] : 

S lF [q,q'}=-2 dt / dt' A(t)V(t,t')Q(t') + - / dt / dt' A(t)Af(t, t')A(t'), (B7) 
Jti Jti 2 J t . J t . 

where A(t) = q(t) - q'(t) and Q(t) = [q(t) + q'(t)}/2. 
The kernels can be computed as 

V{t, i') = ^ J gz(p) Tr [p [fc ( _p)(t), <MO]) ' ( B8a ) 
AA(M') = y / ^Xb)Tr(p{^ I( _ p) (t)^ Ip (i')}), (B8b) 

where (f>\{x, t) is the field operator in the interaction picture and <fii p (t) is the p-mode of the 
same field operator. By integration by parts, the influence action can also be expressed as 

Si F [q,q'}= P dt 1' dt' A(t)H(t,t')Q(t')+ l - f * dt t dt! A(t)N(t, t')A(t'), (B9) 
Jti Jti 2 J t . J t . 

or as 



S 1F [q,q'}=-2 ['dt [ dt' A(t)D(t, t')Q(t') + / ' dt 5n 2 A(t)Q(t) 
Jti Jti Ju 

+ % - £ dt J* dt' A(t)N(t, t')A(t'), 
where the different kernels are defined as 



(BIO) 



H(t,t') = -2^[8(t-t')V(t,t')} (Blla) 

= -26{t - t')D{t - t') + 5& 2 5(t - t'), (Bllb) 

D(t,t') = ~V(t,t'), (Bile) 

N(t,t') = ^M(t,t'). (Blld) 
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The kernels D(t,t') and N(t,t') are called respectively dissipation and noise kernels. The 
frequency shift 5ti 2 is a formally divergent quantity given by 



xo2 or dV(t,t') 

5SI = 2 hm 



t-*t' 



Ot 



(Bile) 



The dissipation and noise kernels can be computed following Eqs. (B8) and (Bll). The 
value of the dissipation kernel in the frequency space is [68] : 



D H = -s-zV). 



(B12) 



The dissipation kernel is closely related to the kernel H(u) [see Eqs. (Blla) and (Bllb)], 
which is also state-independent and given by: 



y ' y J 2vr uj - uj' + ie 



(B13) 



In contrast, the noise kernel (Blld) is state-dependent. For a general Gaussian stationary 
environments, characterised by the occupation numbers n(p) = TV (p e apO p ), the noise kernel 
in Fourier space is given by 



N(u>) 



g 2 \uj\2(u) 



+ n M 



The dissipation and noise kernels are related through: 

~1 



N(u) = — isign(a;) 



+ n(M) 



D(u). 



(B14) 



(B15) 



For the particular case of an environment in thermal equilibrium at a temperature T the 
occupation numbers are given by = l/(e^' T — 1) and the above equation becomes 

the fluctuation-dissipation theorem: 



N(u) = —i sign(w) coth 



\L0\ 

2T 



D(u>). 



(B16) 



By considering an arbitrary distribution of frequencies and an arbitrary Gaussian 
state for the environment p e the dissipation and noise kernels may adopt almost any value. 
In the rest of the appendix we shall try to express all results in terms of the dissipation and 
noise kernels. To this end, it will prove useful to reexpress Eq. (Bllb) in Fourier space: 

dJ iD(u' 



+ 8W 



2n uj — uj' + ie 

The kernel H{uj) can be decomposed in its real and imaginary parts as: 



H r (uj) =ReH(uj) = -2P 
Hi{uj) =lmH(uj) =iD(uj 



duj' iD(uj') x ~ 2 
2tt uj — uj 



(B17a) 

(B17b) 
(B17c) 
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We have used the property l/(x + ie) = P(l/x) — m8{x). Notice that H(—lo) = H*(u>) = 
H R (u) — iH\{u)). The real and imaginary parts of the kernel in frequency space correspond, 
respectively, to the even and odd parts in time space. Notice also the Kramers-Kronig 
relation between the real and imaginary parts of the kernel H{uj): 

AJ Hi(uj') 



-2P 



+ 5& 2 . 



(B18) 



2ir uj — uj 1 

The frequency shift can be always absorbed in O. From now on, and in the main body of 
the paper, we will assume that such absorption has been carried out. 

For a Gaussian environment and asymptotic initial boundary conditions the generating 
functional can be expressed as 5 [83] 

Z[j ll j 2 ] = e^ d * ldt2di3d * 4jA(tl)GR( * 1 -* 2)iV(t2 '' 3) - ?A( * 4)GR( * 4 ' i3) e ~/ dild * 2jA ( tl ) GR (* 1 '* 2 )' ?E (* 2 \ 

(B19) 

where = [ji(t) + j2(t)]/2, j'a(^) = ji(t) — j2{t) and G R (t,t') is the retarded propagator 
of the kernel 

d 2 



L(M') = 
i.e., is the kernel which verifies 

dsG R (t,s)L(s,t'] 



dt 2 



+ n 2 ^j 5{t-t') + H{t,t'), 



-iS{t-t'), G R (t,t') = iit<t'. 



(B20) 



(B21) 



Explicit expressions are most easily obtained in the Fourier space, in which the retarded 
propagator adopts the form: 



— i 



(B22) 



L(w) -uj 2 + n 2 + H(u) ' 

It can be checked that Gn(t, tf), besides being the retarded propagator of the kernel L(t, t'), 
in the sense of Eq. (B21 ), is also the retarded propagator of the quantum mechanical system, 
in the sense of Eq. (Alia). 



Differentiating the CTP generating functional, according to Eqs. (A8), we obtain the 



different correlation functions in terms of the noise and dissipation kernels: 

-i [-uj 2 + n 2 + H r (uj)] + N(u) 



G F (u) 
G + {u) 



[-u J 2 + n 2 + H R (u J )} 2 + [H 1 (u;)} 2, 
i [-uj 2 + 9? + H k (oj)) + N(w) 
[-u 2 + £l 2 + H R {u;)] 2 + [H l {u)] 2) 

N{u) + Hjju) 

[-u; 2 + n 2 + H R (u;)} 2 + [H l (tu)} 2, 

N{uj) - Hj(uj) 
[-cv 2 + n 2 + H R (uj)} 2 + [H l (u;)] 2 ' 



(B23a) 
(B23b) 
(B23c) 
(B23d) 



In the general case there would be a prefactor taking into account the initial conditions for the system. 
However, since the system has a dissipative dynamics and the initial conditions are given in the remote 
past, initial conditions for the system turn out to be completely irrelevant. 
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The generating functional can be alternatively expressed in terms of these correlation func- 
tions as 

Z[ji,j 2 ] =Ne-^ dtdt '^ G -»^ b ^ . (B24) 

If desired, the dynamics of the Brownian oscillator can be analyzed in terms of a Langevin 
equation, 

/•oo 

q(t) + tfq(t) + / dt' H(t, t')q(t') = f (i), (B25) 
Jti 

where £(i) is a stochastic Gaussian field defined by the correlation functions 

<£(i)> 5 = 0, (tmt')h = W,t% (B26) 

with (■■■)£ meaning stochastic average. The stochastic correlation functions derived from 
the Langevin equation correspond to a subset of the quantum correlation functions [3] [83] . 

It is also possible to study the master equation for the reduced density matrix of the 
system p s , which is given by [60J [831 EZJ: 



2 \dq 2 dq' 2 J 2 



+ r(t)/(^-,0(| + | 7 



p s (q,q',t) 



where Stl^t), T(t), h(t) and f(t) are a frequency shift, a dissipative factor and two dispersive 
factors, respectively, which in the weak coupling limit are given by [60] : 



5n 2 (t) = -2 J D(s, ti) cos 0(s - ti) ds, (B28a) 
1 /•* 

T(t) = - D(a,ti) sin n(a-*i)ds, (B28b) 

r(t)fc(t) = / iV(s,ti)cosO(s-ti)ds, (B28c) 

1 /"* 

r(t)/(t) = - / N(t h s) sin fi(a - k) ds. (B28d) 

The corresponding expressions valid for arbitrary couplings can be found in Refs. [60|, [83] . 
Recall that the time U is the time at which the density matrix is assumed to factorize. 
Alternatively, it is also possible to introduce the reduced Wigner function, 

W B (q,p,t) = ^- [ dAe iA Pp s (q-A/2,q + A/2,t), (B29) 
2vr J 

in terms of which the master equation adopts a Fokker-Planck form: 

dW s ^ s 2 2 9W S d P W s ^ r(Ah(A d 2 W s ^ r(Af(A d 2 W s 

-W = - p ^q- + [n +m ^ )]q ^ +2m ^ +mKt) ^ +mf{t) ^q (B30) 
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The Wigner function has some similarities with a classical distribution function, although 
it cannot be interpreted as a true probability density, among other reasons, beacuse it can 
adopt negative values [3] . 
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